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Abstract 

We present a theoretical model describing magnetar giant flares on the basis of 
solar flare/coronal mass ejection theory. In our model, a preflare activity plays a 
crucial role in driving evaporating flows and supplying baryonic matters into the 
magnetosphere. The loaded baryonic matter, that is called "prominence", is then 
gradually uplifted via crustal cracking with maintaining a quasi-force-free equilibrium 
of the magnetosphere. Finally the prominence is erupted by the magnetic pressure 
force due to the loss of equilibrium triggered by the explosive magnetic reconnection. 
The giant flare should be induced as a final outcome of the prominence eruption 
accompanied by large-scale field reconfigurations. An essential difference between 
solar and magnetar flares is the control process of their evolutionary dynamics. The 
flaring activity on magnetars is mainly controlled by the radiative process unlike the 
solar flare governed by the electron conduction. It is highly suggestive that our model 
is accountable for the physical properties of the extraordinary giant flare observed on 
2004 December 27 from SGR1806-20, including the source of baryonic matters loaded 
in the expanding ejecta observed after the giant burst. 
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1. Introduction 

There has recently been growing evidence that soft gamma-ray repeaters (SGRs) and 
anomalous X-ray pulsars (AXPs) are the same population of the ultra-strongly magnetized 
neutron star (B > 10 14 G), so called "Magnetar" (Duncan & Thompson 1992; Harding & Lai 
2006). Activities in these objects are powered by the dissipation of strong magnetic fields 
unlike rotation-powered pulsars and accretion-powered X-ray binaries. Both SGRs and AXPs 
generally undergo the quiescent phase with persistent X-ray emission and the recurrent phase 
of soft gamma-ray (Mereghetti et al. 2004). Typical luminosities at these two phases are 
L x ~ 10 33 1 35 erg s _1 and L 7 ~ 10 38 — 10 41 erg s -1 respectively (Kaspi et al. 2003; Woods & 
Thompson 2006). 

Besides the common short burst, the giant flare with enormous energy and long bursting 
duration is exceptionally observed from SGRs. These are the most energetic galactic event 
currently known (~ 10 44 -10 erg). Only three of giant flares have been observed: SGR0526-66 
on 1979 March 5 (Mazets et al. 1979), SGR1900+14 on 1998 August 27 (Hurley et al. 1999; 
Kouveliotou et al. 1999; Feroci et al. 2001; Tanaka et al. 2007), and SGR1806-20 on 2004 
December 27 (Hurley et al. 2005; Palmer et al. 2005; Terasawa et al. 2005). 

The giant flare from SGR1806-20 is the most recent and energetic one. It is characterized 
by an ultra-luminous hard spike, with energy ~ 10 46 erg, lasting ~ 0.1 s, which decays rapidly 
into a soft pulsating tail lasting hundreds of second. The spectrum of the hard spike is fitted 
by the blackbody radiation with the temperature ~ 10 9 K. In addition, a preflare activity with 
the total energy ~ 10 41 erg and the duration ~ 1.0 s is detected 142 s before the main burst. It 
also shows the single blackbody spectrum with the temperature ~ 10 8 K (Boggs et al. 2007). 

The exceptional event observed in association with the giant flare from SGR1806-20 
is an expanding radio emitting ejecta (Taylor et al. 2005; Cameron et al. 2005; Gaensler 
et al.2005). The emission properties are well resolved by the synchrotron radiation from the 
shocked baryonic shell with the mass > 10 24,5 g and the expansion velocity ~ 0.4c if it is roughly 
spherical (Gelfand et al. 2005; Granot et al. 2006). However, the origin of the baryon-load 
of the ejecta remains unsettled although it is essential for promoting better understanding of 
flaring activities on magnetars (see Thompson & Duncan 1995, hereafter TD95; Gelfand et al. 
2005; Lyutikov 2006). 

In these situations, the solar flare gives an important prototype context for the astro- 
physical flaring activity. It is well known that solar flares are also accompanied by analogous 
mass ejection events, that is called "coronal mass ejections (CME)". The physical proper- 
ties of the solar flare/CME are very similar to those of the magnetar flare (Lyutikov 2006). 
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Specifically, the initial spike and subsequent tail emissions associated with magnetar flares are 
naturally reminiscent of impulsive and decay phases of solar flares. In this paper, we construct 
a theoretical model for the magnetar giant flare on the basis of the solar flare/CME model. 

Our paper is organized as follows. The magnetic reconnection model, which is the 
underlying theory of the solar flare/CME, is applied to the magnetar system for providing a 
physical basis of our magnetar model in § 2. In § 3, we propose a theoretical model describing 
the magnetar giant flare according to a promising solar flare/CME scenario. In § 4, we discuss 
the effectiveness of the assumptions used in our model and the baryonic evaporation process in 
the microscopic viewpoint. Finally, we summarize the characteristics of our model in § 5. 

2. Magnetic Reconnection Model 

Magnetic reconnection is believed to be a crucial mechanism of the energy release in 
the solar flare (Parker 1963; Petschek 1964). The observational evidence, such as cusp-shaped 
soft X-ray loops and hard X-ray sources above the loops, supports this model that predicts the 
primary site of the energy release above the soft X-ray loops (Tsuneta et al. 1992; Masuda et 
al. 1994). The discovery of escaping plasmoids from the flare sites is another evidence for this 
model because it predicts the ejection of plasma from the reconnection region (Shibata et al. 
1995; Ohyama & Shibata 1997). 

The released magnetic energy triggered by the magnetic reconnection is converted to the 
thermal energy in the magnetosphere and then the heat conduction drives the evaporation of 
chromospheric plasma. The discoveries of blue-shifted component of spectral lines and moving 
plasma in X-ray images, confirm the upward motion of the ablated plasma anticipated from 
the chromospheric evaporation theory (Feldman et al. 1980; Culhane et al. 1992; Doschek et 
al. 1992; Savy 1997). 

Shibata & Yokoyama (1999; hereafter SY99) propose a magnetic reconnection model of 
the solar flare taking account of the chromospheric evaporation. They point out that the flare 
temperature is determined by the counterbalance between the reconnection heating and the 
conductive cooling. The coronal density is then controlled by the evaporation cooling of the 
chromospheric plasma which compensates for the conductive heating (see also, Shimojo et al. 
2001; Miyagoshi & Yokoyama 2003). 

The scaling relation obtained from the magnetic reconnection model can explain the 
observed correlation between the emission measure EM and the flare temperature T from solar 
micro-flares to proto-stellar flares consistently (Feldman et al. 1995; Yokoyama & Shibata 1998, 
2001; Aschwanden et al. 2008). This suggests that the underlying physics of the flare would be 
common to various astrophysical systems. We apply the magnetic reconnection model to the 
mysterious magnetar system for providing the physical basis of the flaring activity. 
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2.1. Energetics of Magnetar Flare 

2.1.1. Energy Release by the Magnetic Reconnection 

We consider a situation in which enormous magnetic energy stored in the magnetosphere 
is released through a flare induced by the magnetic reconnection. Then a lot of magnetic arcade 
loops are formed on magnetar surface. Here the reconnected single flare loop is schematically 
illustrated in Figure la. Figure lb focuses on a reconnection site to clearly specify the energy 
conversion process. 

At first, we provide a physical process that controls the typical energy liberated by 
magnetar flares. When we follow a classical magnetic reconnection model, the energy release 
rate of a flare can be described using the released magnetic energy flux F mag 

c ^ flarc _ p 4 fi\ 

(Priest & Forbes 2002) where E^ arc is the released energy by a flare, and A is the area of the 
reconnection site. 

As shown in Figure lb, the magnetic reconnection can liberate the magnetic energy 
which is equivalent to that inflows into the reconnection site. Given the magnetic energy 
density B 2 /4n and the inflow velocity Vi n , the magnetic energy flux -F mag is 
B 2 

-^mag ^ ^in j (2) 

where B is the strength of the magnetic field. 

If we consider only a single flare loop, the area of reconnection site would be given by 
A = LW, where L is the the height of the reconnection point and W is the width of the single 
flare loop (see fig. la). However, we now suppose that there are formed a lot of magnetic arcade 
loops on the magnetar surface. Hence for taking account of the contribution from all magnetic 
loops, we adopt A ~ LR as the total area of all reconnection sites, where R is the size of the 
active region sustaining flaring activities. We fix the typical size R as 10 6 cm in the following. 

Since the reconnection timescale is evaluated as t rcc ~ L/V{ n , equation (1) becomes 
B 2 

-E'flare = I - -R-^ , (3) 
47T 

where the approximation dEa^/dt ~ E^ c /t TCC is used for obtaining this equation. Based on 
the classical magnetic reconnection theory, the field strength B and the reconnection height L 
mainly controls the magnetic energy liberated by the flare. 

2.1.2. Thermal Balance of Magnetar Flare 

In the case of the solar flare, the flare temperature T is determined by the conductive 
cooling which balances with the reconnection heating (SY99). However, the radiative cooling 
dominates the conductive one in the physical condition realized in the magnetar flare (see § 4.3 
in detail). The temperature in the flaring state of the magnetar is thus mainly controlled by 
the radiative cooling which compensates for the reconnection heating. 



Additionally, we should take account of the heat absorption due to the pair creation 
which would impact on the thermal balance of the magnetar system especially in the higher 
temperature regime T< 10 10 K. Considering that the blackbody cooling and the heat absorption 
due to the pair creation become predominate in various cooling processes, the thermal balancing 
equation gives 



where At is the flare duration, Ob is the Stefan-Boltzmann constant, c is the speed of light, 
m e is the electron mass, and n± is the density of the electron positron pair. Note that U ra d 
represents the radiative energy due to the blackbody cooling and U e ± shows the endothermic 
energy due to the pair creation. The magnetic effects on the thermal equilibrium are expected 
to be vanishingly small and ignored here because the configuration of post-flare loops would 
not be changed by post-reconnection processes. 

We focus on the system with B = 10 15 G and T <C 10 10 K which are suitable conditions 
for describing the flaring activity of the magnetar (Boggs et al. 2007). Assuming the local 
thermal equilibrium (LTE), in the range B 3> -Bqed and T <C m e c 2 /ks — 10 10 K, the density of 
the electron positron pair supplied by pair creations is 



[see eq. (75) in Thompson & Duncan 2001], h = 2nh is the Planck constant, fee is the Boltzmann 
constant and -Bqed = m 2 c 3 /eh = 4.4 x 10 13 G is the magnetic flux density at which the energy 
of the first electron Landau level becomes comparable to the electron rest mass. By solving 
the nonlinear equation (4) coupled with the equation (3) and (5), we can derive the flare 
temperature as a function of four physical parameters, T = T(L,B,R,At) or T(E fla _ Te ,B,R,At). 

2.2. Baryonic Evaporation during Magnetar Flare 

Once the solar flare begins, the released energy is rapidly transported to the top of the 
chromosphere by the electron heat conduction and heats the chromospheric plasma suddenly. 
Then the pressure of the heated plasma increases drastically and drives the upward flow into 
the magnetic loop. A hot post-flare loop, which is filled by the evaporated dense plasma, should 
be formed finally (Hirayama 1974; Sylwester 1996; SY99). It can be anticipated to operate a 
similar baryonic evaporation process in the magnetar flare (see Liu et al. 2002 for an application 
of chromospheric evaporation process to accretion disks). 

In the magnetar flare, the photon flux plays a crucial role in the heat transport unlike the 
solar flare dominated by the electron heat conduction. This is because the mean free path of the 
electron is very short and it thermalizes instantaneously in the magnetar's magnetosphere. Even 
if the energy is transported by the created electron-positron pair, the evaporation eventually 



£ flare = 4nR 2 cAt (U Tad + U e ±) 

= 4:irR 2 cAt aBT 4 /c + m e n±c 2 



(4) 




(5) 
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occurs since the thermal equilibrium should be established by the photon flux after the pair 
beams are thermalized. 

The incident energy flux, which inflows into the crustal surface, should thus counterbal- 
ance with the outgoing enthalpy flux of the evaporation flow (see fig. la). When we consider 
the contribution from the created electron positron pair during the flare, the number density 
of the baryon in the evaporation flow n cv is provided by the balancing equation; 

^heat = {h + h±)v ev , (6) 

(c.f., SY99) where -F heat = _E flare /(47nR 2 At) is the inflowing downward energy flux, v ev is the 
upward velocity of the evaporation flow, h is the specific enthalpy of composite gas of the 
baryon and the equilibrium radiation field, and h± is that of the created electron positron 
pair. Notice that this equation is equivalent to the energy conservation equation for a mass 
conserving steady system with the composite gas and the created electron positron pair (c.f., 
Mihalas & Mihalas 1984). 

Assuming the thermal equilibrium state, the specific enthalpy of the composite gas h is 

h = n ev k B T(-+4a^ , (7) 
a = — — = — , (8) 

-fgas 3C71 

(Mihalas & Mihalas 1984) where n cv is the number density of the baryon loaded in the evap- 
oration flow, and a is the ratio of radiation and gas pressures. The specific enthalpy of the 

electron positron pair h± is additionally given by 

5 

h± = n±m e c 2 + -n±k B T , (9) 

where n± is the pair density given by equation (5). 

The upward velocity of the evaporation flow v cv can be replaced by the sound speed of 
the composite gas C s , according to the chromospheric evaporation theory of the sun (SY99; 
Shimojo et al. 2001), 

1/2 

(10) 



c 



s 



rraevfc jB Tfl are (l -fa) 
n ev m p c 2 + h 



5/2 + 20a + 16a 2 

" (3/2 + 12a)(l + «) ' ( } 

(Mihalas & Mihalas 1984) where m p is the proton mass. We thus find from the equation (6) 
that the baryon density of the evaporation flow is the function of E^ are , B, R, At since the flare 
temperature is given by T = T (Eft aTe ,B,R,At). 

Note that it is not settled whether the significant downward energy flux can be main- 
tained after the radiation and pair energy density sufficiently grows at the surface in the realistic 
magnetar situation. We naively assume here that it is sustained as long as the flaring activity 



6 



lasts. The duration sustaining the strong downward photon flux will be investigated in our 
future numerical work. 

2.3. Numerical Solutions Describing the Evaporation 

For given parameters E aaTC , B, R and At, we can obtain T, v ev , a, n± and n cv by 
iteratively solving the coupled equations (4)— (11). The evaporated baryonic mass M ev is derived 
from a relation M ev = Air R 2 m p n ev v ev At (Shimojo et al. 2001). Note that the size of the active 
region and the field strength are fixed as R = 10 6 cm and B = 10 15 G respectively. 

Figure 2a shows the flare temperature T, number densities of the evaporated baryon 
n cv and the created electron positron pair n± as a function of an arbitrary parameter E^^ in 
the case with the fixed flare duration At = 1.0 sec. The each variables are normalized by their 
typical values. 

As is expected, all physical variables increase with increasing the flare energy. Notice 
that the growth rate of the flare temperature is slightly reduced in the high energy range 
-^flarc ^ 10 44 erg by the cooling effect due to the pair creation. However, the baryon density 
steadily increases because the enthalpy of the created pair is negligible compared with that of 
the baryon in the range of interest. 

The mass of the evaporated baryon M ev , the evaporation velocity v ev , the ratio of ra- 
diation and gas pressures a, and the ratio of energy densities for the electron positron pair 
and the radiation field U e ±/U ISl d are demonstrated as a function of the flare temperature T 
[= T (Ea aTe ,B,R, At)] respectively in Figure 2b. The flare duration is fixed as At = 1.0 sec 
again. 

It is found from this figure that the ratio of energy densities U e ± /C/rad increases with the 
increase of the flare temperature, and reaches to unity around T ~ 10 9 K. This indicates that 
the pair production cooling partially contributes to the thermal balance of the system in the 
range T < 10 10 K. 

Furthermore, the ratio of the radiation and gas pressures a slightly decreases in the high 
temperature range. This is because the increasing rate of the radiation pressure is suppressed 
due to the pair creation although the gas pressure increases steadily when T > 10 8 K. The 
apparent enhancement of the growth rate for the baryonic mass M ev reflects the reduction of 
the temperature increasing rate due to the pair production cooling. 

We finally depict the evaporated baryonic mass M ev in Figure 3 as a function of the flare 
temperature T in the cases with different flare durations At = 0.1, 1.0, and 10.0 sec. We find 
that the evaporated baryonic mass is definitely larger in the case with longer flare duration. 
The characteristics of curves are almost same in three models. Numerical results indicate that 
the pair production effect does not dramatically change the qualitative features of our model 
in the temperature range T <C 10 10 K. 

We can thus neglect the contributions from the radiation field and the created pair 
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plasma in the evaporation process described by equation (6). Since the thermal balance of the 
system can be also retained only by the blackbody cooling in equation (4), we can analytically 
derive the flare temperature and the evaporated baryonic mass using scaling relations, 

T ~ (lQ7r 2 Ra B Aty l/A B^ 2 L^ 2 , (12) 

and 

M ev * ^^aTk^m p R^E^ e At^ , 

- ^a B kz l m p R 2 T 3 At . (13) 



These can be used only when the flare energy -E^rc <C 10 48 erg and the flare temperature 
T < 10 10 K. 

3. Solar-type Magnetic Reconnection Model for Magnetar Giant Flare 

A promising solar flare/CME scenario, which is strongly supported by observational 
and theoretical studies (e.g., Priest & Forbes 2002; Shibata 2005), indicates that the baryonic 
material loaded in the CME is evaporated from the sub-coronal chromospheric region before 
its erupting stage. The evaporated matter is then trapped in the coronal region and held in the 
mechanical equilibrium retained by the balance between the magnetic tension and magnetic 
pressure forces. The CME event is finally driven by the magnetic pressure force after the 
loss of equilibrium which is caused by the dissipation of the magnetic tension via the magnetic 
reconnection. The CME induces large-scale field reconfigurations and triggers the main bursting 
activity as the final outcome. 

According to the solar flare/CME scenario, we propose a model describing magnetar 
giant flares on the basis of the underlying magnetic reconnection theory constructed in § 2. 
Our model consists of the following four stages which are illustrated in Figure 4 schematically. 

(a) . A flaring activity begins from preflare stage. In this stage, the magnetic energy 
stored in the magnetar's magnetosphere is partially liberated by the magnetic reconnection 
(Fig 4a). Supposing that the typical height of the reconnection point is relatively low and is 
given by an order of L = 10 3 cm, the released energy and the temperature characterizing the 
preflare, E pre and T pre , are evaluated from equations (3) and (12), 

E PIC = 8.0 x 10 40 B 2 15 R 6 L 2 3 erg, (14) 
T pre = 1.0 x lO 8 ^ 2 ^ 1/4 Lf At - 1/4 K , (15) 

where B^, Ato, Rq and L3 are the field strength, the flare duration, the size of the active 
region, and the reconnection height in units of 10 15 G, 1.0 sec, 10 6 cm, and 10 3 cm respectively 
(c.f., Boggs et al. 2007). 

(b) . During the preflare activity, the radiative heat flux transports the released energy 
and heats the crustal sub-surface matter. The pressure of the crustal matter then increases 
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drastically and drives the upward evaporation flow. As the result, a hot and dense prominence, 
which is trapped by the post-flare loops, is builded up (fig. 4b). The mass of the prominence 
M pro is comparable to that of baryonic matters evaporated and is given by equation (13) 

M pro = 3.4 x 10 24 i? 2 T 8 3 At g, (16) 

where T 8 is the normalized preflare temperature T/10 8 K [see the typical preflare temperature 
in eq. (15)]. 

Immediately after the preflare stage, the formed prominence is bound to the lower mag- 
netosphere of the height ~ 0(100) cm gravitationally, while the reconnection point where the 
preflare is triggered is an order of 10 3 cm [see eq.(14)]. This is because the potential energy 
of the evaporated matter becomes restricted by the input radiative energy (liberated magnetic 
energy) . 

(c) . After the preflare stage, the prominence is gradually lifted up by the magnetic 
energy injected from the magnetar's interior during quiescent stage (see § 4.1 for details of the 
prominence uplifting process). The system evolves with retaining a quasi- force- free equilibrium 
by counterbalancing the magnetic tension with the magnetic pressure. The prominence is 
finally erupted by the magnetic pressure force after the loss of equilibrium which is caused by 
the dissipation of the magnetic tension via the magnetic reconnection. The prominence eruption 
induces large-scale field reconfigurations and triggers a giant burst as the final outcome (Fig 4c). 
Supposing the reconnection height for liberating numerous magnetic energy at the main burst 
stage to be L ~ 4 x 10 5 cm, the energy and temperature of the giant burst -E ma in and T main are, 
from equations (3) and (12) 

£ main = 1.3xl0 4 ^B 2 5 i? 6 L 2 erg, (17) 
T main = 3.7 x 10 9 £} 5 /2 i? 6 - 1/4 Lf At} 74 K , (18) 

where L 5 = L/(4 x 10 5 cm) and At_i = At/10^ 1 s (c.f., Hurley et al. 2005; Terasawa et al. 
2005). 

(d) . The released energy at the main burst stage should be converted into the kinetic 
energy of the erupted prominence and the radiative energy of the remained flare loops. The 
ejected baryon-rich prominence, accelerated by the main burst, would be observed as an ex- 
panding ejecta (Fig 4d). We would like to emphasize again that the baryonic matter loaded 
in the ejecta is supplied by the evaporation at the preflare stage, that is M C j = M pro ~ 10 24 g. 
This is consistent with the observed value of the baryon load in the radio emitting ejecta in 
association with the giant flare from SGR1806-20 (c.f., Gelfand et al. 2005; Granot et al. 2005). 

The remained flare loops are, in contrast, polluted again by the baryonic matter evapo- 
rated by the main burst. The mass of the evaporated baryon in this stage M main is also given 
by equation (13) 

M main = 1.7 x 10 28 R 6 TiAt^ g . (19) 
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where Tg is the temperature at the main burst stage normalized by 3.7 x 10 9 K [see the reference 
value of eq. (18)]. 

As will be discussed in §4.2, the evaporated baryonic matter at the main bust stage 
would be gravit at tonally trapped to the magnetar surface. The baryon-rich dense flare loops 
would be the origin of a trapped fireball, and which eventually produce a luminous 7-ray spike 
and subsequent hard X-ray pulsating tails as is observed in the giant flare from SGR 1806-20. 

Our solar-type magnetic reconnection model is accountable for the flaring activity asso- 
ciated with the giant flare on 2004 December 27 from SGR 1806-20 consistently. An important 
prediction from our model is that the preflare activity plays a crucial role in supplying the 
baryonic matter into the potential ejecta "prominence". This suggests that the radio afterglow 
is expected to be observed only after the giant flare with preflare activities, such as the giant 
burst of SGR 1806-20. 

We would like to stress that the mechanism for baryonic eruptions proposed in our 
model is magnetic pressure-driven one which is caused by the loss of equilibrium triggered by 
the magnetic reconnection. The preflare- induced mass evaporation plays a role in supplying the 
baryonic matter into the magnetospheric region. The mass ejecting mechanism in our model 
is thus essentially different from the magnetic tension-driven model via slingshot like process 
proposed by Gelfand et al. (2005). 

4. Discussion 

4-1. A Possible Process for Uplifting the Prominence 

We discuss a possible physical mechanism for uplifting the prominence during the quies- 
cent stage in the magnetar system. In our model, we suppose that the magnetic energy required 
for uplifting the prominence is supplied from the magnetar's interior via the crustal cracking 
by the Lorentz force like as the model proposed by Lyutikov (2006). This is because the dipole 
field of the magnetar is strong enough to deform the neutron star crust. 

During the quiescent stage, the crustal surface would be deformed by the Lorentz force 
and the magnetic energy stored in the magnetar's interior is converted to the motional energy 
of the crust. The crustal deformation induces the twisting of magnetic fields attached to the 
magnetar surface and generates helical field components. The magnetic energy and flux injected 
into the magnetosphere then re-configure the magnetospheric field. The evaporated baryonic 
matter (= prominence), which is gravitationally bound to the magnetar surface just after the 
preflare stage, would be lifted upwardly in association with the field reconfiguration in the 
magnetosphere . 

According to Lyutikov (2006), when we consider a crustal plate of size R rotating under 
the influence of the Lorentz force, balanced by viscous stress at the base of the curst, the 
dissipated energy by the crustal cracking can be evaluated as ~ 10 44 i? 4 / o 14 (T rot /0.1 sec) -1,5 [erg], 
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where pu is the density normalized by 10 g cm -3 and T rot is the rotation period of the deep 
crustal plate (Landau & Lifshitz 1975). Note that the dissipation energy depends on the 
rotation period of the crustal plate T rot . 

Although the relation between the typical deformation time T rot of the deep crust and 
the duration of the quiescent stage is not settled clearly yet, the crustal cracking by the Lorentz 
force should be a promising mechanism, which is alternative to the buoyant flux emergence in 
the case of the sun, for supplying the magnetic energy into the magnetosphere. 

The process for lifting the prominence remains largely speculative. It is our future 
work to clarify the physical process for triggering the prominence eruption in the magnetar 
system. For verifying the validity of our model, we are now working on the systematic study 
of magnetar's flaring activities using relativistic MHD simulations (Matsumoto et al. 2010 
submitted). 

4-2. Suitable Stage for the Baryon Loading 

We discuss the suitable stage for supplying the baryonic matter into the potential ejecta. 
There are two candidate stages, one is the preflare stage and the other is the main burst stage. 
From equations (16) and (19), the baryonic masses evaporated during each stages are 

M prc ~ 3.4 x 10 24 T 3 # 2 At g , (20) 
M main ~ 1.7 x 10 28 T|i2jjA*_i g , (21) 

where M pre is the evaporated mass during the preflare stage which is comparable to M pro . The 
observational constraint on baryonic mass loaded in the ejecta from SGR 1806-20 is M ej > 10 245 
g (Gelfand et al. 2005; Granot et al. 2006). Both stages can supply the sufficient baryonic 
matter satisfying the constraint. 

On the other hand, the gravitational binding energies of the evaporated matters are, at 
each stages, 

E g , pT e = GM^sM pre /R^s 

~ 6.8 x 10 44 T 3 i2 2 At erg , (22) 

-Eg,mam = G MNsM main / i?NS 

~ 3.4 x 10 48 T 3 i2 2 At_ 1 erg , (23) 

where E g ^ re is the binding energy of the baryonic matter supplied during the preflare stage, 
-Eg.main is that during the main burst stage, G is the gravity constant, M NS = 1.5M is the mass 
of the neutron star, and i?NS is the radius of the neutron star given by 10 6 cm. 

It is found that the binding energy of the evaporated matter supplied during the main 
burst stage is much larger than the bursting energy of the giant flare given by equation (17), 
that is -Emain — 10 46 erg <C -E^main- The baryon supplied by the main burst should be trapped 
on the magnetar surface without escaping. On the other hand, the preflare- supplied baryonic 
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matter can escape, by accelerating the giant burst, from the gravitational field of the magnetar. 
It should be thus the preflare activity that supplies the baryonic matter into the potential ejecta. 

4-3. The Process Sustaining the Thermal Equilibrium 

4-3.1. Optically Thick Flare Loop 

In our model, we naively assume that the blackbody cooling mainly retains the thermal 
balance of the system. Here we validate the effectiveness of this assumption. At first, the optical 
thickness of post-flare loops should be examined to check the availability of the blackbody. 

We can obtain the number densities of the evaporated baryon during the preflare and 
main burst stages n pre and n main from equations (20) and (21) 

- M pre /{m p LR 2 ) 

= 2.1 x 10 33 L3 x 7f At cm" 3 , (24) 

^rnain — -^main/ {jThpLR ) 

= 1.1 x 10 34 L6 1 r 9 3 At_ 1 cm" 3 . (25) 

The Rosseland mean scattering cross-section in the direction parallel to the magnetic field is 

a cs = 2.2 x 10 9 T 2 B~ 2 a T , (26) 

(Silantev & Yakovlev 1980) where o~t is the Thomson scattering cross-section defined by 
(87r/3)(e 2 /m e c 2 ) 2 . Using equations (24)-(26), the optical depths of post-flare loops are given, 
at each stages, 

Tpre TlpreO~esl- J 

~ 3.1 x 10 7 T|.Bf B 2 Afo , (27) 

''"main ^main^es-^ 

~ 2.1 x lQ u T*B^At-i . (28) 

These indicate that the post-flare loops are optically dense and can be treated as the blackbody 
sources. 

4-3.2. Dominant Cooling Process 

Using physical parameters describing the preflare stage, we compare energy evacuation 
rates in various cooling processes and confirm that the blackbody cooling plays a main role 
in retaining the thermal balance of the system. The cooling rate sustained by the blackbody 
radiation from optically thick post-flare loops is 

Abb = o- B T A 

~ 5.7 x 10 27 T 8 4 erg cm _2 sec _1 . (29) 

There are two other cooling processes expected in the magnetar's magnetosphere. One 
is the radiative heat conduction A r and the other is the electron heat conduction A e (TD95). 
The cooling rates by each conductive processes are 
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A r = K r dT/dz 

~ 1.4 x lO^Tg-^Lg 1 At erg cm" 2 sec -1 , (30) 
A e = K e dT/dz 

~ 1.0 x lO^Tj^Lg 1 erg cm" 2 sec" 1 , (31) 



where 



K r = 16a B T 3 /(3Y e na cs ) , k £ = k T 3/2 . (32) 

Here kq ~ 10 -6 cgs is Spitzer's thermal conductivity and the approximation dT/dz ~ T/L is 
used for deriving equations (30) and (31). These clearly indicate that the blackbody cooling 
becomes predominant in the preflare stage. Our model can be thus applicable to this stage. 
Note that, at the main burst stage, the blackbody cooling also mainly retains the thermal 
equilibrium of the system. 

4-4- Microscopic Model of Baryonic Evaporation 

We finally discuss the baryonic evaporation from the microscopic view point. In our 
model, the baryonic material is supposed to be heated by collisions between incident high- 
energy photons and the crustal matter. Considering that the ions are located in the center of 
Wigner-Seitz cell of the crust, the total internal energy per nuclei is naively given by the sum 
of Coulomb lattice energy ei at and thermal energy £th 

e tot — £ lat + £ th 

=1.6 x l^Zl(-\% 
+1.4 x 10~ 10 T CTj6 erg/nuclei , (33) 

(e.g., Shapiro & Teukolsky 1983) where Z26 = Z/2Q is the atomic number normalized by that 
of iron, n cri 28 = ^cr/10 28 cm -3 is the normalized crustal density, and T CTj6 = T cr /10 6 K is the 
normalized crustal temperature. 

Since the lattice energy becomes predominant in the dense crustal surface, the total 
number of the evaporating nuclei N should satisfy the following equation; 

A 

-Eflare ^ A^lat ^ / 4:7lR 2 n CT (.Z^lat dz , (34) 
JO 

where I is the traveling depth of the incident photon through the magnetar crust, and z is the 
depth from the magnetar surface. 

The hydrostatic balance between the surface gravity gn = g/10 14 cm s -2 and the 
pressure gradient force of a degenerate relativistic Fermi gas in the strong magnetic field 
P e ~ 7r 2 h 2 c 2 nl ust /{3eB) (TD95, Appendix A) gives 

dP e 

= p(z)g ~m p n cimt (z)g , (35) 
and leads to the density distribution of the crust, using the approximation dP e /dz ~ P e /z, 
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n cmst (z) ~ 2.4 x 10 31 fii 5 ^i4 z 3 , (36) 

where 2:3 is the depth from the magnetar surface normalized by 10 3 cm. 

Substituting the equation (35) into (34), we can obtain the traveling depth of the photon 
through the crust; 

/ = 4.6 x 10 3 E 3 a il A1 R^ /7 B^ /7 g^ /7 cm . (37) 

The amount of the surface material heated by the collision between the high-energy photon 
and the crustal matter is, therefore, 

rl 

M cv = 4irR 2 m B / n crust (z) dz (38) 

J 

_ UY in24 c.2/7 D-l/7 -1/7 p 6/7 

- 5.4 x 1U K 6 B 15 g u i% arCi41 g . 

This is consistent with the value of equation (16) which is derived from the macroscopic view 
point in the framework of the solar-type magnetic reconnection model. Microscopic evaporation 
model would also support our magnetar flare model constructed in § 2 and 3. 

Here we assume naively that all the released radiative energy is spent only for exciting 
ions, not for increasing their thermal and potential energies. Actually, it is not the easy task 
to clearly specify how the released radiative energy is distributed to each energy components 
during the complicated flaring activity. If we follow the energy equipartition law, a fraction of 
the released energy can be spent for exciting ions at least, and is enough for evaporating the 
baryonic matter with the mass O(10 24 ) g. In order to draw a physical picture for the energy 
distribution process more precisely, we need further study on the nature of strongly magnetized 
crystal (Harding & Lai 2006). 

We would like to stress that the magnetic field stronger than 10 14 G can modify structure 
of crystal of ions (e.g. Harding & Lai 2006; Hansel et al. 2007). The binding energy of crystal 
in strong magnetic fields becomes different from that in weak fields. However, the structure of 
crystal in strong fields has not yet been studied well (see however, Usov et al. 1980). The main 
purpose of this paper is to construct a theoretical model for the magnetar's flaring activity 
from the macroscopic view point on the basis of the solar flare/CME model, not to establish a 
microscopic basis for our model. Thus, in this study, we give rough order-estimation of binding 
energy of crystal using a simple "Wigner-Seits" model in equation (33). 

5. Summary 

According to the magnetic reconnection model which can correctly capture the solar 
flare/coronal mass ejection event, we propose a theoretical model for magnetar giant flares. It 
is highly suggestive that our model is accountable for the flaring activity associated with the 
giant flare from SGR 1806-20 consistently. Our main findings and characteristics of our model 
are summarized as follows: 

1. The temperature of the magnetar flare is essentially determined by the radiative 
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cooling which compensates for the reconnection heating. The cooling effect due to the pair 
production partially contributes the thermal balance of the system, but is not significant. Since 
the blackbody cooling retains the thermal balance of the system, the flare temperature can be 
given by a simple scaling relation T oc i? -1 / 4 At -1 / 4 B l l 2 L 1 / 2 . 

2. During the flaring activity, the photon flux transports the released energy and heats 
the magnetar crust. Then the pressure of the crustal matter increases drastically and drives the 
evaporation flow. The incident radiative heat flux balances with the outgoing enthalpy flux of 
the evaporation flow. Neglecting the enthalpy contribution from the electron positron pair and 
the radiation field, the mass of the evaporated baryon can be represented by a scaling relation 
M ev cx R 2 T 3 At. 

3. In our model, the preflare activity plays a role in supplying the baryonic matter into 
the magnetar's magnetosphere. The "prominence" which contains preflare- supplied baryonic 
matters is gradually uplifted via the energy injection from the magnetar's interior with main- 
taining a quasi-force-free equilibrium of the magnetosphere. Finally the prominence is erupted 
by the magnetic pressure force due to the loss of the equilibrium triggered by the magnetic 
reconnection at the main burst stage. The giant flare should be induced as the final outcome 
by the prominence eruption accompanied by large-scale field reconfigurations. 

4. Our model predicts that the preflare activity produces a baryon-rich prominence. 
Then the erupted prominence is the origin of the observed radio-emitting ejecta associated 
with the giant flare from SGR 1806-20. In contrast, the post-flare loop formed in the main 
burst stage is polluted by dense baryonic matters and be trapped to the magnetar surface. 
This should be the origin of trapped fireball which causes the ultra-luminas 7-ray spike and the 
pulsating X-ray tail in the extraordinary flare observed from SGR 1806-20. 
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Fig. 1. Panel (a): A single magnetic arcade loop formed after the magnetic reconnection which triggers 
the explosive magnetic energy release. We suppose now the situation in which a lot of magnetic arcade 
loops arc formed on the magnctar surface. The typical reconnection height and width of a single flare loop 
is represented by L and W. The radiative heat flux driven via the magnetic reconnection heats the crustal 
matter and drives the upward evaporation flow into the flare loop. Panel (b): A schematic view which 
focuses on the reconnection site. The energy release rate by the magnetic reconnection is comparable to the 
inflow rate of the magnetic energy into the reconnection site, that is Ea aIe /t Iec — F mag A = (B 2 /4tt) V\ a A. 
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Fig. 2. (a): The flare temperature T, the evaporated baryon density n ov and the created electron positron 
pair density n± as a function of an arbitrary parameter Eft aie in the case with the fixed flare duration 
At = 1.0 s. (b): The mass of the evaporated baryon M ev , the evaporation velocity v cv , the ratio of the 
radiation and gas pressures a and the ratio of energy densities for the electron positron pair and the 
radiation field U e ±/U la d as a function of the flare temperature T [ = T(Efi aic , B, R, At)] in the case with 
At = 1.0 s. The size of the active region and the field strength are fixed as R = 10 6 cm and B = 10 15 G in 
these two charts. Note that all physical variables are normalized by their typical values. 
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Fig. 3. The mass of the evaporated baryon M ev as a function of the flare temperature T in the cases 
with different flare durations At = 0.1, 1.0 and 10.0 sec. The size of the active region and the field strength 
are fixed as R = 10 6 cm and B = 10 15 G here. The mass of the evaporated baryon becomes larger when 
the longer flare duration. 
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Fig. 4. A theoretical model for magnetar giant flares on the basis of the solar-type magnetic reconnection 
theory: (a). At preflarc stage, enormous magnetic energies stored in the magnctosphere is partially 
liberated by the magnetic reconnection. The radiative heat flux transports the liberated energy and heats 
the magnetar's crust. The pressure of crustal matters then increases suddenly and drives evaporating flows, 
(b). A baryon-rich prominence with a mass M pro ~ 10 24 g is builded up due to the mass evaporation in 
the lower part of the magnetosphere [~ O(10 2 ) cm], (c). The prominence is gradually uplifted by the 
magnetic energy supplied from beneath the magnetar surface via crustal cracking during quiescent stage. 
Finally, the prominence eruption, which is initiated by the loss of equilibrium triggered by the magnetic 
reconnection, induces large-scale field reconfigurations and triggers a giant burst, (d). The liberated 
energy at the main burst stage is converted to the kinetic energy of the erupted prominence and the 
radiative energy of trapped fireball. The ejected prominence would be observed as radio emitting ejecta 
and the trapped fireball is expected to be the origin of 7-ray spike followed by pulsating hard X-ray tails. 
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